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^ ' Abstract: 

Tj- ! The associated production of a J/ip and a tt meson in antiproton-nucleon annihilation is 

studied in the framework of QCD colhnear factorization. In this approach, a hard subpro- 
cess responsible for the production of the heavy quark-antiquark pair factorizes from soft 
hadronic matrix elements, such as the antiproton (nucleon) distribution amplitude and 

kJ • the nucleon-to-pion (antiproton-to-pion) transition distribution amplitude. This reaction 

H ' mechanism should dominate the forward and backward kinematical regions, where the 

" " " cross sections are expected to be measurable in the set-up of the PANDA experiment at 

the GSI-Fair facility. 
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1 Introduction 

The factorization of exclusive amplitudes in a short distance dominated part which may 
be calculated in a perturbative way in the one hand, and universal confinement related 
hadronic matrix elements in the other hand, is a welcome feature of quantum chromo- 
dynamics (QCD) for specific processes in specific kinematics. The textbook example 
of such a factorization is nearly forward deeply virtual Compton scattering where gen- 
eralized parton distributions (GPDs) are the relevant hadronic matrix elements. The 
extension of this description to other processes such as backward virtual Compton scat- 
tering and backward meson electroproduction, has been advocated [|l], |^, |^ - although 
not proven. In the latter process, new hadronic matrix elements of three quark operators 
on the light cone, the nucleon-to-meson transition distribution amplitudes (TDAs), ap- 
pear which shed a new light on the nucleon structure, and enable to quantify nucleon's 
mesonic cloud. The validity of such a factorization requires at least the existence of a 
large scale Q, which has been taken as the spacelike, respectively timelike, virtuality 
of the photon quantifying the electromagnetic probe in the case of electroproduction, 
respectively lepton pair emission in antinucleon nucleon annihilation. This large scale 
ensures the validity of the perturbative expansion of the hard subprocess with the QCD 
coupling constant OgifJ') taken at a scale fi = 0{Q). The factorization scale is also to be 
taken of order Q. 

In many instances (jet production, heavy meson production or decays...) it has been 
shown that the occurrence of a heavy elementary particle, for instance a quark with mass 
rriQ, is sufficient to ensure the reliability of a perturbative expansion with /i = mg. The 
pioneering studies of the charmonium decay width and of specific charmonium decay 
channels 0, §|, §] have proven the value of this approach. We will follow this line of 
reasoning in the present paper. 

The spectroscopy of charmonium states is at the heart of the physics program of 
the PANDA experiment at the GSI-Fair facility 0, ^. Some of these unknown (or badly 
known) states will decay in a lower mass charmonium such as the J/ijj and a few ordinary 
mesons. When scrutinizing the final states produced, it will be of the utmost importance 
to separate these associated charmonium-light meson states from a background where the 
light meson(s) just evaporate from the beam or the target nucleon. These latter cases 
are the target of our study. But we stress that their production process is much more 
interesting than a mere background to subtract: it is indeed a new way to access the inner 
structure of the nucleon through the study of baryon to meson transition distribution 
amplitudes. For definiteness, we will present the study in the case of a single vr meson 
produced in conjunction with the J/ip. A slight extension of the same formalism applies 
in the cases where another meson, e.g. p, t], oj, /q, if ot a. pair of mesons, e.g. vrvr, KK, is 
produced. 

An alternative description of NN annihilation in a charmonium accompanied by a 
meson proposed in the literature is the use of effective hadron exchange models; see 
e.g. and |10|, |l^. This method provides an essentially non-perturbative description 
of the process. Our aim is to single out with the help of factorization techniques the 
perturbative part of the process and to relate the remaining non-perturbative part to 
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universal fundamental quantities with interpretation within QCD as the light-cone matrix 
elements of correlators of the fundamental fields. 



2 Kinematics 

The study of charmonium exclusive hadronic decays has been for a long time one of the 
fields of application of perturbative QCD methods. It has been argued [|l^, || that the 
dominating mechanism is the cc pair annihilation into the minimal possible number of 
gluons which then produce quark-antiquark pairs forming the outgoing hadrons. 

In the present paper we extend the same perturbative QCD framework for the de- 
scription of a cross channel reaction in which nucleon-antinucleon annihilate producing 
heavy quarkonium together with a light meson (vr, 77, w, p) almost collinear either with 
incoming nucleon or antinucleon. For simplicity below we address the case of vr meson 
and consider the reaction 

N{j>^) +N{p^) -^J/^{p^) + rt{p^). (1) 

Here the NN center-of-mass energy squared s = {pn +Piv)^ — ^^ cind the charmonium 
mass squared Ml introduce the natural hard scale. In the complete analogy with our 
analysis |T3|, |14| we assume that this reaction admits a factorized description within two 
distinct kinematical regimes (see Fig. |l]): 

• the near forward kinematics t = {p.,^ —PnY ~ Oj i^ corresponds to the pion moving 
almost in the direction of the initial antinucleon in NN center-of-mass system 
(CMS); 

• the near backward kinematics u = {Ptt—PnY ~ corresponding to the pion moving 
almost in the direction of the initial nucleon in NN CMS. 

The suggested reaction mechanism forms the distinct forward and backward peaks of the 
differential cross section of the reaction (|lD da/dt as the function of cos 6'* (6'* is the pion 
scattering angle in the NN CMS). The process (|l]) allows us to test the universality of 
T^N TDAs that appear also in the description of 7*A^ — )■ t^N and NN — )■ t^l'ii reactions 

mm. 



Due to the C-invariance of strong interaction there exists a perfect symmetry between 
the forward and backward kinematics regimes of the reaction (||) . These two regimes can 
be considered in exactly the same way (see the discussion in Appendix C of Ref. [Il4| ). 
Precisely, the amplitude of the reaction (|l]) within the t-channel factorization regime can 
be obtained from that within the w-channel factorization regime with the obvious change 
of the kinematical variables: 

Pn -^ Pn, Pn ^ Pn] 

A" = (Pn - Pn) ^ a* = {p^ - pfj); 

u^t. (2) 
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Figure 1: Collinear factorization of the annihilation process N(pff)N{pjsr) — )• J /ip{p^)'K{p-,^). 
Left panel: backward kinematics (u ~ 0). Right panel: forward kinematics (t ~ 0). N{N) 
DA stands for the distribution amplitude of antinucleon (nucleon); itN{ttN) TDA stands for 
the transition distribution amplitude from a nucleon (antinucleon) to a pion. 



However, the backward and the forward regimes are treated somewhat unequally 
within the PANDA experimental set-up operating antiproton beam. Indeed, once switch- 
ing to the laboratory system (which corresponds to the nucleon at rest) one may check 
that the forward peak of the cross section as the function of cos 6'^^^ is narrowed, while 
the backward peak is broadened by the effect of the Lorentz boost from the NN CMS 
to the laboratory frame. 

In this paper we have chosen to present explicitly the details of calculation of the 
reaction amplitude within the backward kinematics regime. As usual, the z axis is chosen 
along the colliding nucleon-antinucleon with the positive direction defined by that of the 
antinucleon beam. We introduce the light-cone vectors p,n satisfying 2p ■ n = 1. The 



Sudakov decomposition of the relevant momenta is presented in the Appendix A of 1 14 



For the calculation of the hard part of the amplitude we apply the collinear approxi- 
mation. We neglect both the nucleon and pion masses and assume A^ = {pn —Pn)t = 0, 
where the transverse direction is defined with respect to the z direction. 

The Sudakov decomposition employed for the calculation of the hard part then reads: 



Pn 



[1 + Op; 



p^ ~ 2^p + 



Pn 



-n; 
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Here ^ is the w-channel skewness variable 

{Pn - Pn) ■ n 



;i+o 



-2^p. 



(3) 
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{p. + PN)-n 2W^-Ml 



Following [12 , in our calculation we set 
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M, 



(4) 



(5) 



taking the average value M = 3 GeV. 



Let us emphasize that we keep the exact kinematics (see App. A of Ref. |T^) for the 
nucleon spinors and for J /if) polarization vector. The physical kinematical domain for 
the reaction (0) in the backward regime is determined by the requirement A|. < 0, where 
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where ttitv {f^-n) stands for the nucleon (pion) mass. The limiting value A\ = corre- 
sponds to 



A^ 



At.. = 



In the calculations presented below we neglect the pion mass 171^^. 



(7) 



3 Hard amplitude calculation for N + N ^ J/iJj + tt 

The calculation oi N + N ^ J/ip + vr scattering amplitude follows the same main steps 
as the classical calculation of the J/ip -^ p + p amplitude |jl2|, |, ^ . Assuming the 
factorization of small and large distance dynamics the hard part of the amplitude is 
computed within perturbative QCD. Large distance dynamics is encoded within the ma- 
trix elements of QCD light-cone operators between the appropriate hadronic states. The 
leading order amplitude of (|T]) is then given by the sum of the three diagrams presented 
on Figure 0. 
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Figure 2: Feynman diagrams describing J/?/' vr production subprocess at the Born order. 



Below we present our conventions for the relevant light-cone matrix elements encod- 
ing the soft dynamics. For definiteness we consider the case of the leading twist uud pn^ 
TDA. The case of nn~ TDA is completely analogous. Throughout this paper we make 
use of the parametrization of |[T7[| P| since it simplifies considerably for A^ = case. The 
parametrization involves 8 invariant functions each being the function of three longitudi- 
nal momentum fractions Xi, skewness variable ^, momentum transfer squared A^ as well 



^The relation between this parametrization and the one employed in Refs. 11^, IT^, Q is given in the 
Appendix A of [^. 



as of the factorization scale /i: 
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+m^2Tf^")(xi,2,3,e,A2)(a,A,C),.(A^t/+), 

Here /^ = 93 MeV is the pion weak decay constant and f^ determines the value of 
the nucleon wave function at the origin. Throughout this paper we adopt Dirac's "hat" 
notation v = v^'y^] a^" = ^[y^,'y'^]', cr^^ = Vxa^^; C is the charge conjugation matrix and 
f/+ = pn U{pn, Sn) is the large component of the nucleon spinor. For A^- = just three 
invariant amplitudes V^ , A^^ and T}^^ survive in the parametrization (H). 

For the leading twist antinucleon DAs we employ the standard parametrization [^ 
(see also Appendix B of Ref. [0)- The non-relativistic light-cone wave function of J/ip 
heavy quarkonium is given by p[ 



'^pr{z,p^) = {0\Cr{z)Cp{-z)\J/llj) = - f^ 



2mcS + (Tp^^S" 



(9) 
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where rUc is the c-quark mass and S stands for the charmonium polarization vector. With 
the use of the non-relativistic wave function (j^) we tacitly assume that each charm quark 
carries half of the momentum of the J/ip. The normalization constant /^ is extracted 
from the charmonium leptonic decay width T{J /ip — )■ e^e~): 
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Using the values quoted in [^ we get 

l/^l =413 ±8 MeV. 
The leading order amplitude of (P reads 
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where V and U stand for the nucleon Dirac spinors. The calculation of 3 diagrams 
presented on Fig. yields the following result for X(^, A^): 

J-^+^ i=l -^0 fc=l 

\ ymy^ixi + ^0)(X2 + z0)(x3 + z0)(xi(2yi - 1) - 2^^! + zO) (0:3(2^3 - 1) - 2^1/3 + ^0) 






+ 



yiy2y-i{xi + 20)(x2 + i0)(x3 + i0)(xi(2|/i - 1) - 2^1/1 + i0)(x2(2y2 - 1) - 2^1/2 + iO) 

(13) 
where 

/ d^x= i dxi dx2 dx3] dsy = / dyi dy2 dy^. (14) 

J-1+5 -/-i+C ^-i+C -/-i+c Jo Jo Jo Jo 

ForX'(^,A2) we get 



X'(e,A2)= / rf3X(5(Vx,-20 fd^ydiS^yu-l 



i=i " '^ fc=i 

[ e^(Xil/3 + X3yi)(r2^'^^Hxi,2,3,e,A^) -4-^^X1.2,3, e, A^))(\/ni/l,2,3) -^ni/1.2,3)) 

\ ymy^ixi + «o)(x2 + 20)(x3 + 20)(xi(2yi - 1) - 2^1/1 + 20)(x3(2|/3 - 1) - 2^1/3 + ^o) 

e^(Xll/2 + X2yi)(rf^Hxi,2.3,e,A^)+rf^)(Xi,2,3,e,A^))Tnyi,2,3) 

yi2/22/3(a;i + z0)(x2 + ^0)(X3 + «0)(xi(2yi - 1) - 2^2/1 + ^0)(x2(2y2 - 1) - 2^2/2 + ^0) 

(15) 
The overall factor C in (0) is expressed as: 

n tA \3JNJ^ 1 ,, ic X. //I \3JNJi' 10 /ic\ 

C = 47ra» — -— - X 16 x -■-■—-—-= 47raj — -— — , 16 

^ '^ /^ 2 ,'^^ 3 3 (3!)2 ^ '^ f^ 81' ^ ^ 

'"^,7^' '. color factor 

normalization 

where as stands for the strong coupling. One may check that the poles in Xi are indeed 
located either on the cross over trajectories Xi = 0, which separate the DGLAP-like and 
ERBL-like support regions of tcN TDAsQ, or within the DGLAP-like support region (for 
the y-dependent poles in Xj) as it certainly should be. 

The structure of result ([T^ , (|13D , (|15D resembles much the well known expression for 
J/ip — !■ pp decay amplitude 



M = {4r.asr^^^USVMo, (17) 



^For the definition of the ERBL-hke and DGLAP-hke support regions of TDAs see 
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(19) 



The J/tp — )■ pp decay amphtude (|T7D resuhs in the following expression for the decay 
width |1|: 

r(JM -> pp) = (TraJ^l^^-^lMop. (20) 

4 Estimates of the cross section 

The squared amplitude (12) averaged over spins of initial particles reads 

\Mxx'\' = \J2 ^T'HMim*- (21) 

SpfSfj 

At the leading twist only the transverse polarization of J/ip is relevant. To sum over the 
transverse polarization we employ the relation: 

5^£^(A)£*'^(A) = -g^'' + ^—{p^n'^+p'^nn, (22) 



and get 



lA^P - E \M^\' = \n^^ (|X(e, A^)P - ^l^'(^'^')l') • (23) 



The leading twist differential cross section oi N + N -^ J/^p + tt then reads ^1 

'dE? " 167rA2(s,m2„m^)'-^^' ' ^^^^ 



where k{x,y,z) = \J x^ + y"^ + z^ — 2xy — 2xz — 2yz. 

In order to get a rough estimate of the cross section we use the simple nucleon exchange 



model for ttN TDAs suggested in ||15[. We do not expect that the inclusion of the spectral 
part for vrA^ TDAs ||2^ would be essential to draw a conclusion on the feasibility of the 
relevant experiment and may be postponed until the precise experimental data will be 
available. 



For the vrA^ TDAs within the parainetrization (§) the nucleoli pole model of Ref. [|T^ 
reads 
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Here V^, A^ and T^ stand for the nucleon DAs; Qt^nn 
nomenological coupling and 
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ensures the pure ERBL-like support of TDAs. For the simple nucleon pole model 
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where Mq is given by eq. (P^. 

As the phenomenological input for our cross section calculation we may use different 
solutions for the leading twist nucleon DA. Similarly to the case of charmonium decay 
width our result depends strongly on the form of input nucleon DA. In order to be able 
to characterize this dependence we present in Table |I| the predictions of J/ijj decay width 
{^^ into pp within the pQCD description with the use of the nucleon DAs in question 
as the numerical input. Several points are to be mentioned. 

• The decay width ( pUD shows strong dependence on Og: ~ a^. There is no unique 
opinion in the literature on the value of the strong coupling for the gluon virtuality 
in question. 
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for as = as = 0.3 


is reproduced 


COZ 


(5.0 ±0.5) -10-3 


0.42 


0.26 


KS 


(5.0±0.5)-10-=' 


0.84 


0.24 


"Heterotic" model 


(5.0±0.5)-10^=' 


1.20 


0.22 


BK 


6.64-10-3 


0.05 


0.38 


BLW NLO 


(5.0 ±0.5) -10-3 


0.02 


0.44 


Asymptotic 


(5.0 ±0.5) -10-3 


0.015 


0.46 


Experiment 


- 


0.19^0.21 


- 



Table 1: Nucleon light-cone wave function normalization constant /at at /io = 1 GeV^ 
and J/ip — > pp decay width for various nucleon DA models. Experimental value is taken 
from [111. 



Nucleon DAs that are strongly concentrated in the end-point regions such as Chernyak- 
Ogloblin-Zhitnitsky (COZ) [|16] or King and Sachrajda (KS) ||2^ for as = 0.3 seem 
to overestimate the experimental width by the factor 2 -=- 4. They require smaller 
values of as ~ 0.25 to reproduce the experimental value. These solutions have been 



strongly criticized in the literature (see e.g. discussion in Chapter 4 of Ref. p3 



The "heterotic" DA model |Q requires even smaller values of as 
the experimental width. 



0.2 to reproduce 



On the other hand, the phenomenological solutions for the nucleon DA which are 
close to the asymptotic form (Bolz-Kroll (BK) [|25|, Braun-Lenz-Wittmann (BLW 
NLO) model of [^) for as = 0.3 underestimate the experimental value of the decay 
width by a factor 5-i-lO. Usually, a^ ~ 0.4 is required to reproduce the experimental 
value of T{J/ip — )• pp). 



Based on the above arguments we have chosen to present our results for the NN — )• 
J/ipTT cross section with the value of a^ fixed by the requirement that the given phe- 
nomenological solution reproduces the experimental J/ip — )■ NN decay width. On Fig. ^ 
we show our estimates of the differential cross section -^ for pp — )■ J/ip 7r° as a function 
of W^ for A'^ = 0. On Fig. § we show the differential cross section ^^ for pp — t- J/ip 7r° 
as a function of A^ for several values of W^. 

The pion scattering angle in the NN CMS for the w-channel factorization regime then 
can be expressed as: 
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Figure 3: Differential cross section -^ for pp — )■ J/tJj n^ as a function of W"^ for A|. = 0. 
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Figure 4: Differential cross section j^ for pp — )■ J /i}) 7r° as a function of A|> for W"^ = 15 
GeV^ (left panel) and W^ = 20 GeV^ (right panel). 
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4(1 + 

One may check that for Ay = indeed cos 6'* = —1, which means backward scattering. 
Expressing our cross sections as functions of cos ^* (^Uj) we reproduce the order of mag- 
nitude (~ 100 -=- 300 pb/GeV^ for A"^ = 0) of the corresponding cross sections presented 
in |ll|. This fact is certainly not surprising since it can be seen as an artefact of a simple 
nucleon pole model of nN TDAs. 

We also reproduce the characteristic shape of angular distributions of the cross sec- 
tions presented in [^ . On Fig. ^ we show the center of mass angular distribution for the 
da/dA'^ cross section for both forward and backward factorization regimes presented on 
the polar plot with the polar angle being the pion CMS scattering angle 6**. We present 
the ratio 



lS^iW',Al) 



(31) 



as the function of 6* showing the result for W'^ = 15 GeV^ and for —1 GeV^ < A^ < A^ 



10 



where A^^^ > is the hmiting value (^ of the momentum transfer squared. The left 
half of the graph corresponds to the near-backward factorization regime and right half 
of the graph corresponds to the near- forward factorization (see Fig|I|). With the dashed 
lines we show the effect of the cutoff A^ = — 1 GeV^ for the values of the CMS scattering 
angle. 

Since these rates are certainly within the experimental reach of the PANDA exper- 
iment, the study of reaction (|1|) will provide a valuable universality test for the TDA 
approach since the same TDAs also arise in the description of NN — > 
backward pion electroproduction off nucleon 7*A^ — )• vrA^ |l^ . 
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It is worth mentioning that it might be advantageous also to study the process 
P{Pn) + Hpn) -^ J/i^iPi,) + 7r"(p^). 



(32) 



In our simple nucleon pole model for t^N TDAs the corresponding cross sections are 
enhanced by the factor 2 due to the isotopic factor \/2 in (EUf) . 
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Figure 5: Angular distribution for the da jdS} cross section for near forward (cos0* > 0) 
and near backward (cos^* < 0) scattering regimes for —1 GeV^ < A^ < A^^^^. Dashed 
lines show the effect of the cutoff A^ > — 1 GeV^ for the values of the pion CMS scattering 
angle Q^. 



5 Conclusions 



In this paper we address the reaction p + A^ — > J/^ + vr which will be studied in the 



PANDA experiment at GSI-FAIR to look for exotic charmonium states production |27 
We argue that outside the region specific for the resonance production, this reaction 
may be analyzed within the pQCD framework.lt will not only help to quantitatively 
disentangle resonance production from the universal hadronic background but also will 
provide valuable information on hadronic structure encoded in nucleon-to-pion TDAs. 
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Nucleon-to-pion TDAs are essentially non-diagonal matrix elements of QCD light cone 
operators which probe the non-minimal Fock state contents of hadrons. Therefore TDAs 
supply complementary information with respect to diagonal partonic distributions (PDFs, 
GPDs). On the other hand the suggested possibility of description of the process in terms 
of fundamental degrees of freedom of QCD largely increases its theoretical importance. 
It is also worth mentioning the possible generalization of our approach both to the case 
of other heavy quarkonium states as well as to various accompanying light meson species 
(?7, p, cj, etc.). The TDA framework has also been recently p8[ used in a double handbag 



description of proton-antiproton annihilation into a heavy meson pair. 

Within the kinematical range accessible at PANDA we provide the predictions using 
a simple nucleon pole model for nN TDAs. The obtained values of cross sections give 
hope of experimental accessibility of the reaction. Our predictions are consistent with the 
recent estimates of []1T| obtained within a fully non-perturbative effective hadronic theory. 



However, the latter approach lacks the direct relation to the dynamics of fundamental 
degrees of freedom of QCD. Precise experimental data, and the study of polarization 
observables not discussed here, will allow to discriminate between the QCD and hadronic 
approaches. 

It is also worth mentioning that the mass of the charm quark may not be large 
enough for our leading order (in a^) and leading twist analysis to be sufficient to describe 
the data. More work is certainly needed to go beyond the Born approximation for the 
hard amplitude, in particular because the timelike nature of the hard probe is often 



accompanied by large 0{as) corrections p9 
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